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Abstract Second-order two-scale expansions, a unified proof for the regularity of the correctors 
based on the translation invariant and a lemma for extracting O(e) from the remainder term arc 
presented for the second order nonlinear elliptic equation with rapidly oscillating coefficients. If the 
data are smooth enough, the error of the zero-order (energy) in first-order in the Holder norm, 
second-order's gradient (flux) in the maximum norm(linear periodic case), are locally O(e). It can be 
. used in the parabolic equation. 
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> • 1 Introduction 

■ ' 

Consider the homogenization of the following elliptic problem: find u E € Hq(Q), 

d ( x du \ x 

-—^(u £ , X ,-)^J=f(u e , X ,-), inil; (1) 

^ ! 

where fl C M" is a bounded Lipschitz domain and the summation convention is used. A £ = (oy ) 
• is symmetric and positive definite; aij(u e , x, y) are 1-periodic in y. In the combined conduction- 
radiation heat transfer! 1 !, = kij(x, |) + Au^bij. Assume all of the data are smooth enough. 

It is difficult to solve this problem numerically because of the rapidly oscillating coefficients. 
One method is the two-scale homogenization. The error estimate in L°° (0(e)) was presented 
by J. L. Lions ct al ^ and LkJ 3 !. Oleinik et al^ proved the 0{e 1 / 2 ) estimate in H 1 . Su et 
all 5 l investigated the quasi-periodic problems; Zhang and CW 1 ' gave a numerical example for 
the Rosseland equation. There are also some other famous methods, such as Multiscale Finite 
Element Mcthod(MFEM[ 6 l) and Heterogeneous Multiscale Method(HMMl 7 ]). 



OO 
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2 Second-order two-scale expansions 

The periodic cell Y = [0, 1]™; W^ er (Y) C {ip £ H\Y) : J y cp = 0} consists of the functions 
with the same traces on the opposite faces of Y. As the same as the linear case^l, we look for 
a formal asymptotic expansion of the form 

x x 

u e (x) = u (x) + eu 1 (x, '-) + e 2 u 2 {x, -) + ... (2) 

£ £ 

where ui(-,y), U2(-,y) is Y-periodic in y. Let y = | , then -J^ — > + Assume the 

coefficient A £ in (1) has the form A £ = A(uq,x, -) + £Ai(uq, U\, x, -) + £ 2 A 2 + ... where A = 
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{aij{uQ, x, |)) is symmetric positive definite and A\,Ai are uniform bounded; the righthand 
side has the form f(u £ ,x, |) = f(uo,x, f ) + e/i + ... with f\ G L°°. Substituting (2) into (1) 
and equating the power-like terms of e, we introduce the following auxiliary functions to make 
the term of order e _1 equal zero 

aij(u ,x,y) 9 ^ m = - f a mj (u ,x,y)^-, V<p(y) G Wp er (Y),l < m < n; (3) 

where N m (uo,x,y) G Wp er (Y)(uo, x are parameters). Then Hi = NidiUo- The problem for the 
part of order e° admits a unique solution iff there exists u Q G Hq(^) sucn * na * 

= J^f(uo,x,y)dy, a?(uo,x) = JjfHj(v. Q ,x,y) + au(u ,x,y)^^-]dy. (4) 

This system is well-posed because of the compensated compactness! 8 !. V<p(y) G W^ er (F), 

9 , ,9 



Ojj(uo, £, y)j—M k i(u ,x, y)j—<p(y) 
oyi oyj 

dN t f dN t , 



^(Z/) - / UkmNt- — . (5) 



'''' 7^7^" = -] Y {ail ^ + Al ' ik + Al ^ ) dy- i 



+ 



y <9Xi 



dN k f dN k ] 
aik + a-u^, / (aik+au— — ) ip{y). (6) 

oyi Jy oyi J 

f(uo,x,y)- f(uo,x,y)]tp(y). (7) 



If <p(y) = 1, the righthand sides of the above equations equal zero. So M k i, Qk,R & Wp er (Y)(uo, 
x are parameters). Let 1*2 = Mkid^ua + QkdkUo + R to make the O(e ) term equal zero. 

3 Error estimate in L°° and C a 

For the regularity of the following periodic problem(the abstract Eq. from(3),(5)-(7)), it can 
be unified by the translation invariant. HY z = Y + z,zE R", find M G Wp er (Y z ) such that 

A(y)VM ■ = J B(y) ■ + d(y)<p, G W* er (Y z ), (8) 

where A, B, d is Y-periodic and J Y d = 0. Then M = N in Wp er (y) after a periodic extension 
if N is the solution to (8) in the cell Y. So the boundary estimate for N is equivalent to the 
interior estimate for M. If B G L q , d G i« /2 , q > n, then N G C Q (Th 8.24^1). 

Notice that Lemma 1.6^ is right for the case v = 1, u G W 1 ^, p > 1. If g(x,y) G L(f2 x 
K"),/ y g = 0, then by the help of W^" 1 ^ and Meyers estimate(Th4.2l 2 ]) there exists 

~$ e {x) G LP(n;R n ), s. t. g(x, -) = div^ E , H^U^R") < Ce. (9) 

Let Z £ = u e - u= u £ -u ~ eui - e 2 u 2 . [A(u E ) - A(u)]Vu = Z e A'(£)\7u = Z £ b(x, §). Then 
d ( , x.dZ E , . x. r7 \ d9i , 
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because u £l u are known(Thl0.7[ 9 l). Z £ = -eu x - e 2 u 2 is 0(e) in L°°(dn) and 0(e 1_a ) 
in C Q (dft). If IMUINUII^II^IMI^ < C,q > n, then ll^lU^Jc^n,) < Ce by 
the maximum principle and De Giorgi-Nash estimate(Th8.16, Th8.24t 9 l), Q' CC f2. Since 
||«i||ooj ||«2||oo < C 5 IK - "olloo < Ce. Consequently, for the energy, \f n A e Vu e ■ Vu e - 
f n A°X7uo ■ Vuo| ^ Ce. If |e«2 He 3 ^ C, then \\u e — uq — euiHcWQ') ^ Ce. For the case of 
nonlinear parabolic equation with mixed boundary conditions, sec I 10 !. 

If in (10) ciij = aij(f), ay G C 7 (F)(or picccwisc smooth), 6 4 = 0,||-,7? € L ,l+,5 (0),(5 > 0; 
Q' is a open set, f2' C f2. Then by the help of Lemma 16' 3 ', sup|y|VZ e | ^ Ce. If |V K M2| ^ 
Ce -1 , then sup^ | V(u e — uq — £Wi)| ^ Ce. It's also true for the tensor case. For the flux, 
supjp- 1 A E V{u e — Uq — £Ul)\ ^ Ce. 

Let w e (x) = [Z £ (ex) — Z e (0)]/e^ 3 \ then consider the following equation 

-div[A(ex,x)Vw E (x) + ew E (x)b(ex,x)} = eF(ex) + Z s (0)div[b(ex,x)]. (11) 

where F(£) = effK£) + ei l{C)- Differentiate both sides formally, by De Giorgi-Nash estimate 
[dw e } c -,(yjT^ < Ce, since |Z e (0)| < Ce. Then locally [V(u e — u — sui)]ci ^ Ce 1-7 . 

For the gradient estimate in parabolic case, see Li and Lit 11 !. Their work was based on the 
piecewise smooth coefficients so it's very important for both theory and practice. 
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